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Abstract 

We calculate electromagnetic neutral pion production off tri-nucleon bound states ( 3 H, 3 He) at 
threshold to complete one-loop order in the framework of chiral nuclear effective field theory. We 
show that the fourth order two-nucleon corrections to the S-wave multipoles at threshold are very 
small. This implies that a precise measurement of the S-wave cross section for neutral pion produc- 
tion off 3 He allows for a stringent test of the chiral perturbation theory prediction for the S-wave 
electric multipole Eq + u . 



#1 email 
* 2 email 
#3 email 
* 4 email 



lenkewitz@hiskp.uni-bonn.de 
evgeny.epelbaum@rub.de 
hammer@hiskp.uni-bonn.de 
meissner@hiskp.uni-bonn.de 



1 Introduction 



In the absence of free neutron targets, light nuclei like the deuteron or three-nucleon bound states like 
3 H (triton) or 3 He can be used to unravel the properties of neutrons. For a recent review on extracting 
the neutron structure from electron or photon scattering off light nuclei, see Ref. pQ. Of particular 
interest in this respect is threshold neutral pion photo- and electroproduction off the nucleon. This 
is one of the finest reactions to test the chiral dynamics of QCD, see Ref. [2] for a recent review. 
Arguably most striking is the counterintuitive chiral perturbation theory prediction (CHPT) that the 
elementary neutron S-wave multipole Eq_^ is larger in magnitude than the corresponding one of the 

proton, -E'o+ P This prediction was already successfully tested in neutral pion photo- [5] and 

electroproduction off the deuteron [6]. However, given the scarcity and precision of the corresponding 
data, it is mandatory to study also pion production off tri-nucleon bound states, that can be calculated 
nowadays to high precision based on chiral nuclear effective field theory (EFT). This framework 
extends CHPT to nuclear physics (for recent reviews, see [HE])- As it is well-known, 3 He appears to 
be a particularly promising target to extract the information about the neutron amplitude. Its wave 
function is strongly dominated by the principal S-state component which suggests that the spin of 3 He 
is largely driven by the one of the neutron. Consequently, neutral pion photo- and electroproduction 
of tri-nucleon systems ( 3 He and 3 H) was considered in Ref. [9] based on chiral 3N wave functions at 
next-to-leading order in the chiral expansion. Here, we extend this calculation to fourth order in the 
chiral expansion, including consistently all next-to-next-to-leading order contributions. This amounts 
to a complete (i.e. subleading) one-loop calculation. 

Experimentally, neutral pion photoproduction off light nuclei has so far only been studied at 
Saclay pT0]JTT] and at Saskatoon [i~2"l[T3]. Clearly, new measurements using CW beams, modern targets 
and detectors are urgently called for. The results presented below show that we are able to calculate 
neutral pion photoproduction off tri-nucleon systems to very good precision and, moreover, that the 
S-wave cross section for neutral pion production off 3 He is very sensitive to the elementary EqJ 1 
multipole (as already stressed in Ref. [9]). 

Our manuscript is organized as follows. Section 2 contains all the necessary formalism. In particu- 
lar, we spell out in detail the fourth order two-nucleon corrections that modify the third order results 
of Ref. [9]. We also briefly recall the contributions worked out in that paper. Section 3 contains our 
results and the discussion of these. We end with a short summary and outlook in Sec. 4. 

2 Formalism 

2.1 Generalities 

Pion production off a tri-nucleon bound state is given in terms of three different topologies of Feyn- 
man diagrams, see Fig. [TJ While the single-nucleon contribution (a) corresponding to the standard 
impulse approximation features the elementary neutron and proton production amplitudes, the nu- 
clear corrections are given by two-body (b) and three-body (c) terms. Based on the power counting 
developed in [14], at next-to-next-to-leading order (NNLO), only the topologies (a) and (b) contribute. 
Topology (c) starts to contribute at fourth order to P-wave multipoles and is thus not of relevance for 
our considerations. Here, we will specifically consider threshold photo- and electroproduction param- 
eterized in terms of the electric E'en- and longitudinal Lq + S-wave multipoles. In particular, we study 
the sensitivity of the 3 H/ 3 He S-wave multipoles to the elementary -E7g + n multipole, as the production 
amplitude off the proton is well understood experimentally and theoretically. 

To analyze pion photo- and electroproduction of the tri-nucleon system, one has to calculate the 
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Figure 1: Different topologies contributing to pion production off the three- nucleon bound state (tri- 
angle), (a), (b) and (c) represent the single-, two- and three- nucleon contributions, respectively. Solid, 
dashed and wiggly lines denote nucleons, pions and photons, in order. Topology (c) does not contribute 
to the order considered here (NNLO). 



nuclear matrix element of the transition operator O as: 

(M'j\d\Mj)^ := M'j if jv q\6\ij) Mj P 3N k] 



where ip refers to the three-nucleon state and k, q, P37V and P^ N denote the momentum of the real or 
virtual photon, produced pion and the initial and final momenta of the 3N nucleus, respectively. Notice 
that we use the Coulomb gauge throughout this work. The 3N bound state has total nuclear angular 
momentum J = 1/2 with magnetic quantum numbers Mj for the initial and Mj for the final nuclear 
state. J can be decomposed in total spin S = 1/2, 3/2 and total orbital angular momentum L = 0,1, 2. 
The total isospin is a mixture of two components, T = 1/2 and 3/2. While the T = 1/2 component is 
large, the small T = 3/2 component emerges due to isospin breaking and is neglected in our calculation. 
The isospin magnetic quantum numbers are = Mt> = 1/2 for 3 He and Mt = Mjv = — 1/2 for 3 H. 
In this paper, we consider neutral pion production by real or virtual photons off a spin- 1/2 particle - 
either the nucleon or the 3 H and 3 He nuclei. At threshold, the corresponding transition matrix takes 
the form 

M x = 2i E 0+ (?a,t • S) + 2i L 0+ (?a,l • S) , (2) 

with e\,T = £\ — (£\ k)k and e\,L = {c\ k)k the transverse and longitudinal photon polarization vectors, 
and S is the spin vector. The transverse and longitudinal S-wave multipoles are denoted by Eq + and 
Lq + , respectively. Note that Lq + contributes only for virtual photons. Also, as will be explained later 
in more detail, due to boost effects that appear at next-to-next-to-leading order (NNLO), there will 
be P-wave corrections to the S-wave multipoles. In the following, we will first recapitulate the leading 
one- loop calculation from Ref. [9] and then systematically work out the NNLO corrections. 

The chiral expansion of the pertinent amplitudes is done in complete analogy to the case of the 
deuteron discussed in detail in Ref. [5]. We summarize here only briefly the most important issues. 
First, the production operator on the single nucleon is worked out to fourth order in the chiral 
expansion. This is consistent with the two-body contributions (exchange currents) worked out here. 
The leading two-nucleon terms start at third order and have been evaluated in Ref. [9j. At fourth 
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order, there are further corrections to the two-body currents which can be classified as boost, static 
and recoil contributions. These will be worked out in detail in the following paragraphs. Second, using 
the same arguments as in Ref. [5], one can show that there are no contributions from short-distance 
four-nucleon-pion-photon operators at threshold (which is the case we consider here). Notice that 
such operators do contribute at fourth order for nonvanishing momenta of the produced pion. 

2.2 Single nucleon and leading one-loop contributions at threshold 

As explained before, the matrix element Eq. ([1]) receives contributions from one- and two-nucleon 
operators at the order we are working. Consider first the single nucleon contribution, given in terms 
of the 1-body transition operator 1N . After some algebra, one finds 

IN I n/TA . _„•="._. C ( T?*°P T?S+V , j^r°n T?S-V\ , • - . . q ( r j^S+V , r n°n jpS-V 



{M'j\0^\MjU = ie KT ■ S M , Mj (E%>F* +V + E$*F*- V ) + ie x , L • S M , Mj \fi*F$+ v + I^Ff 

(3) 

where F T J^ = F^ L ± Fj>/L and F T ', L denote the corresponding form factors of the 3N bound state, 



3 

F T/L ^*A,T/L • Sm'jMj = ^ ( M j\t\,T/L ' &l\Mj 



2 

F^/l CA.T/L • S m , jMj = | (M'j\e x , T/L ■ a X T{\Mj)^ , (4) 

which parameterize the response of the composite system to the excitation by photons in spin-isospin 
space. Sm'jMj are the corresponding spin transitions matrix elements. In the above equation, {ji) 
denote the spin (isospin) Pauli matrices corresponding to the nucleon %. Furthermore, z refers to the 
isospin quantization axis. 

Using the 3N wave functions from chiral nuclear EFT at the appropriate order, the pertinent 
matrix elements in Eq. @ can be evaluated. Here, we use chiral 3N wave functions obtained from the 



N 2 LO interaction in the Weinberg power counting |15pi6]P 5 l In order to estimate the error from higher 



order corrections, we use wave functions for five different combinations of the cutoff A in the spectral 
function representation of the two-pion exchange and the cutoff A used to regularize the Lippmann- 
Schwinger equation for the two-body T-matrix. The wave functions are taken from Ref. |17pi8j and 
the corresponding cutoff combinations in units of MeV are (A, A) = (450,500), (600,500), (550,600), 
(450,700), (600,700). All five sets describe the binding energies of the 3 He and 3 H nuclei equally well 
(after inclusion of the corresponding three- nucleon force). 

The one-body contribution to the 3N multipoles are given by 



•IN . Km f T n°pjpS+V , TTr n ipS-V 



pIN 

L 0+ ~ ~2~ \ L 0+ L + L 0+ b L ) ■ i 5 ) 

Here, ATiat is the kinematical factor to account for the change in phase space from the IN to the 3N 
system, 

m N + M n m 3N 

^1N = —V7 ~ L092 i 6 

with rriN being the nucleon mass and the mass of the three-nucleon bound state. Throughout, 
we denote the neutral pion mass by M v , but the S-wave multipoles are usually given in terms of 



* 5 The consistency of the Weinberg counting for short-range operators and the non-perturbative renormalization of 
chiral EFT are currently under discussion, see the review [7] for more details. A real alternative to the Weinberg 
approach for practical calculations, however, is not available. 
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nucleus 


3 He 


3 H 


F$ +v 


0.017(13)(3) 


1.493(25) (3) 


Ft V 


1.480(26)(3) 


0.012(13)(3) 


F s+v 


-0.079(14)(8) 


1.487(27)(8) 


F s-v 


1.479(26)(8) 


-0.083(14)(8) 



Table 1: Numerical results for the form factors F T , L . The first error is our estimation of the theoretical 
uncertainty resulting from the truncation of the chiral expansion while the second one is the statistical 
error from the Monte Carlo integration. 



the charged pion mass and also the charged pion mass appears in some of the two-nucleon 

contributions. 

We evaluate the matrix elements for the one-body contribution in Eq. @ numerically with Monte 
Carlo integration using the VEGAS algorithm [19] . The results for the form factors F T ,^ are given in 
Table [TJ The first error represents the theoretical uncertainty estimated from the cutoff variation in 
the wave functions. We take the central value defined by the five different cutoff sets as our prediction 
and estimate the theory error from higher-order corrections from the spread of the calculated values. 
Strictly speaking, this procedure gives a lower bound on the error, but in practice it generates a 
reasonable estimate. The second error is the statistical error from the Monte Carlo evaluation of the 
integrals. It is typically much smaller than the estimated theory error and can be neglected. 

We stress that we follow the nuclear EFT formulation of Lepage, in which the whole effective 
potential is iterated to all orders when solving the Schrodinger equation for the nuclear states. As 
discussed in Ref. |20j . the cutoff should be kept of the order of the breakdown scale or below in order 
to avoid unnatural scaling of the coefficients of higher order terms. Indeed, using larger cutoffs can 
lead to a violation of certain low-energy theorems as demonstrated in Ref. [21] for an exactly solvable 
model. 

The error related to the expansion of the production operator is difficult to estimate given that the 
convergence in the expansion for the single nucleon S-wave multipoles is known to be slow, see Ref. [3] 
for an extended discussion. We therefore give only a rough estimate of this uncertainty. The extractions 
of the proton S-wave photoproduction amplitude based on CHPT using various approximations [22J 
lead to an uncertainty Ai?Q + p ±0.05 x 10 _3 /M 7r +, which is about 5%. The uncertainty of the 
neutron S-wave threshold amplitude is estimated to be the same. Consequently, our estimate of the 
error on the single nucleon amplitude is 5%. 

We now switch to the two-nucleon contribution. In Coulomb gauge, only the two Feynman dia- 
grams shown in Fig. [2] contribute at threshold at third order [5] . Their contribution to the multipoles 
can be written as 

EH = K w (4 a >-lf), 

1$ = K W (F^-F^), (7) 

with the prefactor 

KW = 167r(m 3N +M T )(27r)3F3 « °- 135 fm X 10 ^ ' ® 

The numerical value for i^N was obtained using g A = 1.26 for the axial-vector coupling constant, 
F n = 93 MeV for the pion decay constant, and the neutral pion mass M n = 135 MeV. The transverse 
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Figure 2: Leading two- nucleoli contributions to the nuclear pion production matrix element at thresh- 
old. Solid, dashed and wiggly lines denote nucleons, pions and photons, in order. 



nucleus 


3 He 


3 H 


4 a) - [fm- 1 ] 
if) - Ff [fin-*] 


-29.3(2)(1) 
-22.9(2)(1) 


-29.7(2)(1) 
-23.2(1)(1) 



Table 2: Numerical results for the form factors Fj,L — Fj,j L P arame terizing two-body contributions 
in units of fm . The first error is our estimation of the theoretical uncertainty resulting from the 
truncation of the chiral expansion while the second one is the statistical error from the Monte Carlo 
integration. 



and longitudinal form factors Fj?J L and F^L corresponding to diagrams (a) and (b), respectively, are 



^T/L e A,T/L ' ^M'jMj 



£*A,T/L • (<?1 + <?2)(n • f 2 



T 1 T 2. 



q' 2 



Mj 



and 



P T/L e A,T/L ' ^M'jMj 



3( M' 7 



(n ■ t 2 



T 1 T 2 



1 [(q' - k) ■ (at + a 2 )} [e\,T/L • (9 ' - k/2)] 



[{q'-kY + Ml + ]q 



FT 12 



Mj 



(9) 



(10) 



where q' = P12 — p/2 + k/2 is the momentum of the exchanged pion and p\2 = (^1 — k 2 )/2, p( 2 = 
(k[ — k' 2 )/2 are the initial and final Jacobi momenta of nucleons 1 and 2, respectively. The integral 
for the form factors Frj?) L contains an integrable singularity which can be removed by an appropriate 
variable transformation. Then, the form factors can be evaluated using Monte Carlo integration in 
the same way as the form factors for the single-nucleon contribution. Our results for ijy L — F^/L are 
given in Table El The first error is again the theory error estimated from the cutoff variation in the 
chiral interaction as described above. The second error is the statistical error from the Monte Carlo 
integration which is about half the size of the theory error. 

2.3 Fourth order contributions 



Apart from the fourth order corrections already included in the single nucleon multipoles £"o+ and 
I/0+, there are various fourth order corrections that arise due to the presence of the other two nucleons 
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ai 



a 2 




(as) 




/ s 

i' + > 



a 4 



Figure 3: Subleading static two-nucleon contributions to the nuclear pion production matrix element 
at threshold. The filled circle denotes an insertion from the dimension two effective Lagrangian. For 
further notations, see Fig. [2J 



in the tri-nucleon system considered here. First, there are the boost corrections that arise from the 
observation that within a nucleus the threshold for pion production is shifted compared to the free 
threshold. For light systems as considered here, this essentially induces P-wave contributions as 
detailed in Sec. 12.3.11 Second, there are further corrections to the two-nucleon production operator 
displayed in Figs. [3] and 01 These can be grouped in two categories, namely the so-called static and the 
so-called recoil corrections. The static corrections are shown in Fig. [3l They involve - as the leading 
2N corrections do - static propagators but one insertion from the dimension two chiral effective pion- 
nucleon Lagrangian. The recoil corrections feature corrections to the static propagators with only 
insertions from the leading order (dimension one) chiral Lagrangian. These corrections are most 
conveniently derived in time-ordered perturbation theory. The corresponding diagrams are shown in 
Fig. [H The boxes indicate the regions where the two energy denominators whose 1/mAr-expansion 
generates these corrections can appear (see Ref. [23] for more details) . These two types of corrections 
will be discussed in Sec. 12.3.21 and Sec. 12.3.31 m order. 
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Figure 4: Subleading recoil two-nucleon contributions to the nuclear pion production matrix element 
at threshold in time-ordered perturbation theory. The boxes indicate the regions where the two energy 
denominators to be expanded can appear Further notation as in Fig. El 

2.3.1 Boost corrections to the one-nucleon contributions 

The proton and neutron production amplitudes are calculated in the (N, 7) center-of-mass system 
(cms). The boost of a (3.2V, 7)-cms four- vector p to the (N, 7)-cms four- vector p* has the general form 




7 -7/3 \( P °\ = ( 7(p°-P-p) 



fill 



where P is the projection operator onto the /3-direction, i.e. Ps x = (/3 • x)j3, p\\ = Pg p is the parallel 
part, and p± = (1 — Ps)p the perpendicular part of p = pj| + p± with respect to j3. To determine f3, 
we consider k\ + k. In the (2V, 7)-cms this combination has to vanish, i.e. k* + k* = 0. We have 

kf + k* = 7 + A; ) + P^h + kj) + (1 - P )(h + k) = 0. (12) 

Because linearly independent (even orthogonal) vectors have to vanish separately, i.e. 

(l-P0)(ki + k) = , 
-/3(fc° + k°) + P (h +k) = Q, 

we conclude that 

g _ fci + fc _ fci + g _ # 2 - P3/ 2 + 2 g/3 , s 

V-k° + k °-k ? + q o- ^ _ m _ m? +< + q0 < ) 

Near the static limit we have 



ff_ Pi2-P3/2 + 2g/3 / (p{ 2 -py2-q/3) 2 , 1 , 2 + ^ ] 

= p^-m+m + ( 2) th roS hoi d) ^ 2 -pg/2 + ( 2) 
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and thus 7 = (l — /3 2 ) 1//2 = 1 + O (m^ 2 ). Correspondingly, a general four-vector p M in the (37V, 7)- 
cms transforms to the (TV, 7)-cms as 



P 



" - 7 ( p _ 0. = p _ PX2 ~ Ps/ 2 + 2 ?/3 threshold^ ^0 Pl2 - Pi/ 2 . - 



•p 



-)• p 



p, 



«• , - - - Pl2 -PW 2 + 2 ^/ 3 threshold _ p[ 2 -pz/2 

p = — 7PP + p± + 7P|| = P — 2 P > P — 2 — P 



m N 



m N 



(15) 
(16) 



to first order in the inverse nucleon mass. At threshold, the boosted energy and momentum for the 
photon (A; *, k*) and pion (q *, (f ), as well as the time- and space-like components of the photon 
polarization vector (e°*, e^) are thus given by: 



k°* = k l 



-A - fc A 



(p{ 2 


- m) 


k 




m N 




(Pl2 


~ Pi/2) 


q 




m N 






~ PS/2) 




m N 



k* = k 



m N 



m N 



(P12-P3/2), 



(P12-P3/2) 



3 



eA--^-(Pi2-^/2) = r A . 



(17) 



Since = 0, the polarization vector does not change except for the time component. 

To adopt the results for pion production off nucleons which were calculated in the (N, 7)-cms for 
our calculation on the tri-nucleon, the pion momentum in the (3N, 7)-cms at threshold, q = 0, is 
boosted to the (N, 7)-cms value q* = — {M^ / mN)(p' l2 — pg/2) =: —fi(p[ 2 — P3/2) above threshold. The 
corresponding P-wave contribution reads (using the notation from Ref. |24j ) 

M x = 2i(e x ■ S)(q* ■ k)Px + 2i(e x ■ q*)(k ■ S)P 2 + (e^A ■ [q* X k])P 3 

+ 2i(e x -k)Ck-S)(q*- k)(Pi - P 5 - Px - P 2 ) + 2i(e x ■ k)(q* ■ S)P 5 
=: 2i(e XjT ■ S)(q* ■ k)P x + 2t(? A)T ■ • S)P 2 + (e A , T ■ [q* X k])P 3 

+ 2i(e A , L • S)(q* ■ k)P A + 2*(e A>L • • S)P 5 . 

Close to threshold, the P-wave multipoles Pi behave as Pi « Pj|<f*| = /UPi|p 12 — P3/2I with 



pf = +0.0187 fm 2 , if = +0.0240 fm 2 , if = +0.0013 fm 2 , 
P« = +0.0134 fm 2 , Pf = +0.0234 fm 2 , P£ = +0.0003 fm 2 , 



(18) 



where the numerical values refer to the P-wave low-energy theorems for pion photo- [3J and electro- 
production [25]. Corrections to these theorems are beyond the accuracy of our calculation. 

In analogy to the S-wave case discussed above, we define the corresponding P-wave form factors 



' d M'jMj I F\ 



s/v 



£*A,T • Sm'jMj 



(e*A,T • Sm'jMj 



pS/V 



F. 



s/v 



e*A,L • S_ 



M'jMj 



F 



s/v 



e*A,L • S_ 



■ DM'jMj ) F 5 



s/v 



Mj, 
Mj, 

Mj, 
Mj, 



3 (eV • Sx) ((Pi2 - Pa/2) • k) X s ' v I Mj) . 
3(eA,T • {P'12 - Pi/2)) [k ■ S x ) X s ' v | Mj) , 



■h (^T-m.-^mxk]) x s ' v 



Mj 



3 (e A ,L • Sx) ((p[ 2 - Ps/2) ■ k) X s l y I Mj) , 



3(? A , L -fc) Up[2-pU 2 h-Sx) x s ' v 



Mj 



(19) 
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nucleus 


3 He 


3 H 


pS+v 


+0.004(3)(1) 


+0.339(6)(1) 


pS-V 


+0.338(5)(1) 


+0.002(3)(1) 


pS+V 


-0.015(2)(0) 


-0.011(2)(0) 


Fi-y 


-0.011(2)(0) 


-0.015(2)(0) 


pS+V 
r A 


-0.019(5)(4) 


+0.339(6)(4) 


pS-V 
r A 


+0.337(6)(4) 


-0.021(3)(4) 



Table 3: Numerical results for the boost correction form factors F i in units of [fm~ ]. The first error 
is our estimation of the theoretical uncertainty resulting from the truncation of the chiral expansion 
while the second one is the statistical error from the Monte Carlo integration. 



where the spin and isospin operators refer to nucleon 1. Further, we have introduced the notation 
X s = 1, X v = r-f. The contributions from the other nucleons are accounted for by the overall factor 
of three. 

In terms of these form factors, the P-wave contribution to the 3N-production amplitude takes the 
form: 

3 

B y 008 * « -0.546/i J +1/ Pf + ^ s - v if) , (20) 
i=l 

5 

L^boost ^_ 0M6 ^j2(Ff+ v P[ + Ff- v pA , (21) 

i=A 

where F? ±v = Ff ±fV. These form factors are evaluated using the same Monte Carlo method as 
employed for the S-waves. The numerical values are collected in Tab. Note that F% and F$ come out 
to be consistent with zero and and therefore not listed in the table. As before, the proton contribution 
is dominant in 3 H, whereas the neutron one features prominently in 3 He. 

Notice that in contrast to the single-nucleon corrections, we do not need to employ a special 
treatment for boost corrections to the leading two- nucleon contributions. All 1/mjv-corrections to the 
leading three-body contributions to the production operator needed in the calculations are treated on 
the same footing as described in section 12.3.31 

2.3.2 Fourth order two-nucleon contributions 

Similar to the leading order graphs displayed in Fig. [2j one has additional contributions at NNLO, 
which are displayed as diagram classes 01,02,03,04, and 61 in Fig. [3l These involve the insertion of 
one dimension-two operator from the chiral effective Lagrangian. As in the case of the deuteron, only 
the operators ~ l/2mjv> ~ 9A/^ m N or ~ K p,n contribute, while all contributions proportional to the 
LECs Cj vanish at threshold. All these Feynman diagrams give corrections of the form iTyi = MO\2 
and involve multiples of the generic prefactor M = (egjim^)/ (2F^). They lead to contributions of the 
form 

M=2iKf N (6 l2 + 6 2l )^ (22) 
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with 

< = --— 7 1 „w ^ = ~ K ™ , \, - -0.0104fm 3 ^ (23) 
2N 2 m N 8ir(m 3N + M n ) (27r) 3 2mjv 4mjvM ff M 7r + 

to account for the phase space and normalization. The prefactor of the 2N contributions at order g 4 , 
K 2N is suppressed considerably compared to the corresponding prefactor at order q 3 , K 2 jy. This has 
to be kept in mind when we discuss the numerical results for the fourth order corrections. 

We are now in the position to evaluate the various classes of diagrams, following the calculations 
for the deuteron from Refs. [U[23]. The corrections to diagram (a) in Fig. [2] which are given in 
Fig. E](ai)-(a4) read 

,AW,ai _ n r) 2 \ em N g A e*A,T -vi,-. 



* T 12 =( 1 - 2 fA)^T^3 jv—(n-T2-T{T 2 ){q ■ q + 2q -p l2 ) 



e"A,T • [(q" + 2p/ 2 - • g' + z[g' x fc](l + k v \ 

F 3 7/72" 



■ T NN,a 2 _ em N g A 

U 12 - 75s 7575 in-T2-Tl r 2j) 



e"A,T • f-(g ' + - fc) + ipi x (g ' - k)]) a 2 ■ q' 



■rpNN,a4, 2 em N9A \ ~ -/ - -. f r, A \ 

* T 12 = 9 a— 75| 7572 + M 2 + ( r i • T 2 - rfr|) . (24) 

The corrections to diagram (b) in Fig. [2] which are given in Fig. [3^6i ) read 

*T 12 = - (1 - 2 flA ) ^3 ( - - - m , + ) - - - (n .r 2 -r 1 r 2 )(q -q +2q ■ Pl2 ) . (25) 

with q' as before and q" = q' — k. 

2.3.3 l/m^v-corrections to leading two-nucleon contributions 

Furthermore, there are 1/mTv-corrections to leading two-nucleon diagrams some of which are shown 
in Fig. [2j They have been worked out within the Q-box framework in Ref. [23J. The corresponding 
diagrams a[, b[, b 2 , b' 3 , b' 4 , b' 5 , and b' 6 are displayed in Fig. |H The boxes indicate the regions where the 
two energy denominators whose l/m^r-expansion generates the corrections can appear (see Ref. [23J 
for more details). 

For the expressions of the various diagrams, we use the following abbreviations: 



UJ 



' ~ - lq' 2 + M 2 + and u" = Jq" 2 + M 2 + . (26) 



With that, the correction to diagram (a) in Fig. [2] which are given in Fig. H^a^) read 

NN,a[ em N g A £a,t ; ~ M n ) r , , f _ 

»r u = a/ (a/ + M n ) 2 k ' { ~ 2q ~ 2Pl2 + k){Tl ' T2 " TlT2) ' (2?) 

whereas the various corrections to diagram (6) in Fig. [2] which are given in Fig. H{&i)-(&g) take the 
form 

.rr.NNM em N g A <?i • q"e\,T • (<f ' + <f )( w ' + M n ) _ z 

^ (n-Ti-r^) 

^ ( -q' ■ q' - q' ■ f l2 + 2(q' + p 12 ) ■ k - k 2 q' ■ q' + q' ■ & l2 \ 



Mr w' - w" - M T 
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.rpNN,b' 2 _ em N g A <?i ; q"e\,T ; (g" + 9%^ ~ M*) ,-. _ z z 
12 " 16F3 W V«( W ' + W " + M,)2 {Tl ' T2 TlT2 



x 



2w" {g' • q> + g ' • ft n - 2(q' + # 2 ) • + fc 2 } + (a/ + M w ) { -2(«f + p< 12 ) • fc + P }) , 



.rpNN,b' 3 _ em N g A ai ■ q"ex,T ; (g" + 9%^ ~ M^) ,^ ^ _ * z 
12 ~ 16F| w'w" 3 (u;' + w" + M w ) 2 (w' + a;") 2 2 ^ 

/' {2(g" + • - fc 2 } + 2(u/ + M,) {<f • g' + g' • ^ 12 } 

. T NN,b' 4 = em N g A gi ; g^gyr ; (g" + g%</ + M n ) _ _ z z . 
12 " 16F3 w'a;" 3 (a;' + w" + M w ) 2 (w' + a;'0 2lri ' T2 TlT2j 



x w 



/' {g" • q' + g" • ff 12 - 2(g' +ff 12 ) ■ k + k 2 } - (J - M w ) {q 1 ■ q' + q> • p< 12 } 



.^iViV,^ _ em N g A ax ■ g"e\,T ■ (g" + q')(uj' + ^ _ z 2 

12 " 16^3 W y3( w ' + W "-M,)2 lTl " T2 TlT2 



X 



2w" {g" • q 1 + q' ■ ff l2 - 2{q' + ff 12 ) ■ k + k 2 } + {J - M v ) { -2(g" + # 2 ) • + fc 2 }) , 



.rrNN.b's _ em N g A d\ ■ q"e\^ ■ (g" + q')(u' - M nj ,_. _ 
^ 12 ~^6?| a/^' + M^ (n-ra-nr, 

x (y + w " + M u ){2(g-*' + p' 12 )-£-P}) • (28) 

3 Results and discussion 

We are now in the position to evaluate the nuclear S-wave multipoles. They are given as the sum of 
the one- and two-nucleon contributions given in the previous section, 

E o+ = -Eo+ + -^0+ ) L o+ = L o+ + -^0+ • ( 29 ) 

Combining the leading and subleading corrections to the two-nucleon production operators discussed 
above with the subleading chiral perturbation theory results for the single- nucleon multipoles at 0(p 4 ) 

w 

Ef + P = -1.16 x W- 3 /M n+ , Ell n = +2.13 x 1(T 3 /M^+ , 

iffi = -1.35 x 1(T 3 /M W+ , L f+ = -2-41 x 1(T 3 /M^+ , (30) 

we obtain for the threshold multipoles on 3 He and on 3 H the values listed in Table HI Note that the 
values for the single nucleon multipoles in Eq. (|30p are consistent with the data of [26] and a recent 
calculation in the chiral unitary approach |27j . The extractions of the proton S-wave photoproduction 
multipoles based on CHPT using various approximations show a 5% uncertainty [22]. Consequently, 
we assign a conservative 5% error to the single-nucleon multipoles. 

In Table HI adding the first two columns gives the leading one- loop result from Ref. [9j. The 
fourth order corrections are given separately for the boost of the single nucleon terms (third column) 
as well as the contributions described in sections 12.3.21 and 12.3.31 and referred to as static and recoil, 
respectively. The complete one-loop result can be found in the sixth column. The first error given is 
an estimate of the theory error from higher orders in chiral EFT, the second error is the statistical 
error from the Monte Carlo integration. Notice that the statistical error is negligible compared to 
theory error. The 5% error from the single-nucleon amplitudes discussed above is not included in the 
numbers for the theory error, but appears as the second error of the single nucleon contribution in 
the table. Overall, we find that these fourth order corrections for the electric dipole amplitude Eq + 
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3 He 


IN (g 4 ) 


2N (g 3 ) 


lN-boost 


2N-static (g 4 ) 


2N-recoil (g 4 ) 


total 


E 0+ [W- 3 /M n+ ] 
L 0+ [W- 3 /M w+ ] 


+1.71(4)(9) 
-1.89(4)(9) 


-3.95(3) 
-3.09(2) 


-0.23(1) 
-0.00(0) 


-0.02(0)(1) 
-0.07(1)(1) 


+0.01(2)(1) 
+0.07(7) (0) 


-2.48(11) 
-4.98(12) 


3 H 


IN (g 4 ) 


2N (q 3 ) 


lN-boost 


2N-static (g 4 ) 


2N-recoil (g 4 ) 


total 


E 0+ [W- 3 /M w+ ] 
L 0+ [10- 3 /M, + ] 


-0.93(3)(5) 
-0.99(4)(5) 


-4.01(3) 
-3.13(1) 


-0.35(1) 
-0.02(0) 


-0.02(1)(1) 
-0.07(G)(1) 


+0.01(2)(0) 
+0.07(7)(0) 


-5.28(7) 
-4.14(10) 



Table 4: Numerical results for the 3N multipoles. The first error is our estimation of the theoretical 
uncertainty resulting from the truncation of the chiral expansion while the second one is the statistical 
error from the Monte Carlo integration (with the exception of the IN contributions as explained in 
the text). 



for both tri-nuclear systems come out to be very small, much smaller than in case of the deuteron. 
This can be, in part, traced back to the smaller values of the various form factors (for the boost 
corrections) and also to the small prefactor K% N , cf. Eq. (j23|) . for the two-nucleon contributions. For 
the longitudinal amplitude Lq + , the sum of the fourth order corrections is consistent with zero within 
the uncertainties. This can be understood as follows: First, the boost corrections are proportional to 
the P-wave multipole P4, which is much smaller than the corresponding multipoles Pi, P3 that appear 
in the electric dipole amplitude, cf. Eq. (I18p . Second, there are almost perfect cancellations between 
the static and the recoil contributions for both tri-nucleon systems. These cancellations are accidental 
in the sense that they can not be traced back to any symmetry or small prefactor. In summary, we find 
that the chiral expansion for the S-wave multipoles at threshold converges fast and that the largest 
uncertainty remains in the single nucleon production amplitudes. We also remark that the fourth 
order corrections to the electric dipole amplitudes of the tri-nucleon systems are sizably smaller than 
for the deuteron. 

Now, let us concentrate on photoproduction. The threshold S-wave cross section for pion photo- 
production ao is given by 

I A; I da 



o ail 



= \E 0+ r ■ (31) 

q=0 

In Fig. [5l we illustrate the sensitivity of ao to the single- neutron multipole Eq^J 1 . The shaded band 
indicates the theory error estimated from the cutoff variation as described above and a 5% error in 
E 0+ . As shown above, the uncertainties related to the nuclear effects are of the order of one percent, 
i.e. completely negligible. So our estimate of a 10% uncertainty of the 2N contributions in Ref. [9] 
driven by the analogy to the deuteron case [5] turns out to be much too conservative. The vertical 
dashed line indicates the CHPT prediction £J° n = 2.13 x 10" 3 /M 7r +. Changing this value by ±20% 
leads to changes in ao of about ±30%. Thus, the 3 He nucleus is a very promising target to test the 
CHPT prediction for Eq + h . On the contrary, neutral pion production on 3 H is rather insensitive to 
Eq^J 1 : a variation of in the range ... 3 (in units of 10~ 3 / M n +) changes ao only by 1%. 

Next we compare our predictions with the available data. The consistency of the CHPT prediction 
for the single-neutron multipole with the measured S-wave threshold amplitude on the deuteron from 
Saclay and Saskatoon is well established, see Refs. [2J[5]. The reanalyzed measurement of the S-wave 
amplitude for 3 He at Saclay gives E 3 = (-3.5 ± 0.3) x 10~ 3 /M n + [TOIITT], which is related to a 
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Figure 5: Sensitivity of do for 3 He in units of 10 6 /M^ + to the single-neutron multipole Eq^J 1 in units 
of 10~ 3 /M n +. The vertical dashed line gives the CHPT prediction for £?o+™ an d the vertical dotted 
lines indicate the 5% error in the prediction. The shaded band indicates the theory error estimated 
from the cutoff variation and a 5% error in E^J 3 as described in the text. 



Here, we have approximated the A = 3 body form factor Ta of Argan et al. [TT] by the numerically 
dominant form factor F T ~ for 3 He, cf. Tab. [TJ This results in 



assuming the same sign as for our 3 He prediction in Table 0] In magnitude, the extracted value is 
about 25% above the leading one-loop prediction Eq + = —2.24(11) x 10~ 3 /M 7r +. The discrepancy 
is reduced for the fourth order result Eq + = —2.48(11) x 10 _3 /M 7r + , which is about 12% below the 
experimental value in Eq. (133p . Taking into account the errors, our fourth order result is consistent 
with the experiment of Argan et al. [TT]. Given the model-dependence that is inherent to the analysis 
of Ref. [11] . it is obvious that a more precise measurement using CW beams and modern detectors is 
very much called for. 

4 Summary and outlook 

In this paper, we have presented a calculation of the subleading two-nucleon operators for threshold 
neutral pion photo- and electroproduction off the tri-nucleon systems 3 H and 3 He. This completes 
the one- loop calculation at order 0(q 4 ). The single nucleon contributions to that order and the 
leading third order two-nucleon contributions were already evaluated in the corresponding letter [9]. 
The production operator was evaluated in the framework of chiral nuclear effective field theory, in 



according to 




(32) 



E 0+ = (-2.8 ± 0.2) x 10~ 3 /M, 



(33) 
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line with the earlier calculations for production off the deuteron [6]. To this order, it gets both 
one- and two-body contributions. Here, we have given explicit expressions for the fourth order two- 
nucleon contributions stemming from boost corrections for pion production off a single nucleon (such 
contributions only arise from the fact that in a nucleus the threshold for pion production is lowered 
as compared to the nucleon case), from static 2N contributions with one insertion from the dimension 
two chiral pion-nucleon Lagrangian and from recoil corrections to the pion- and nucleon propagators. 
We used the chiral wave functions of Refs. [T5l[T6j to calculate the S-wave 3N multipoles Eq + and Lq + . 
These wave functions are consistent with the pion production operator. 

We have shown that all corrections at fourth order are very small, a few percent for the tri-nucleon 
electric dipole amplitudes and essentially vanishing for the corresponding longitudinal amplitudes. 
This suppression can be explained by very small boost correction and an accidental cancellation 
between the static and the recoil contibutions. We remark that these corrections are sizably smaller 
than in the deuteron case [5]. 

The theoretical uncertainty associated with the cutoff variation in the employed wave functions 
appears to be small (of the order of 3%). The dominant theoretical error at this order stems from the 
threshold pion production amplitude off the proton and the neutron, which is estimated to be about 
5%. Consequently, we have explored the possibility to extract the elementary neutron multipole EqJ 1 
from a neutral pion photoproduction measurement off 3 He. We found indeed a large sensitivity of the 
Eq + amplitude to Eq_^. Given the very small uncertainty of the nuclear corrections as shown here, 
3 He appears to be a promising target to test the counterintuitive CHPT prediction for Eq^J 1 OH]. 

We have also shown that our prediction for the 3 He S-wave multipole £"o+ is roughly consistent 
with the value deduced from the old Saclay measurement of the threshold cross section [11] . A new 
measurement using modern technology and better methods to deal with few-body dynamics is urgently 
called for. 

There are many natural extensions of this work. They include investigating higher orders, pion 
production above threshold, the extension to virtual photons and pion electroproduction, production of 
charged pions, and considering heavier nuclear targets such as 4 He. Further work in these directions is 
in progress. Last but not least, the role of nucleon recoil, accounting for the new scale \ = yJM^rrifq ~ 
340 MeV, needs to be investigated in more detail in view of the findings of Refs. [28,29j. 
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